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Abstract 

We consider a globally neutral Lorentzian plasma as a possible remnant of a preinfla- 
tionary stage of expansion and pose the problem of the suitable initial conditions for the 
evolution of the large-scale electromagnetic inhomogeneities. During the protoinflationary 
regime the Weyl invariance of the Ohmic current guarantees that the comoving conductivity 
is approximately constant. The subsequent breaking of Weyl invariance by the masses of 
the charge carriers drives the conductivity to zero. The newly derived conducting initial 
conditions for the amplification of large-scale magnetic fields are contrasted with the con- 
ventional vacuum initial conditions. It is shown, in a specific class of examples, that when 
the number of inflationary efolds is close to minimal the effects of the conducting initial 
conditions cannot be neglected. 
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In recent years attention has been paid to the generation of large-scale magnetic fields in 
the early Universe [1]. Of particular interest for the present analysis are those mechanisms 
and toy models whose purpose is to amplify the vacuum fluctuations of the electromagnetic 
fields potentially present at the onset of the inflationary phase or, more accurately, around 
63 efolds prior to the end of inflation. Inflation cannot be however eternal in its past and the 
inflationary stage of expansion is normally complemented by a pre-inflationary stage where 
the evolution of the background is decelerated and, most likely, dominated by radiation 
[2]. Protoinflationary initial conditions different from the vacuum have been discussed by 
various authors [2] in connection with the scalar and tensor modes of the geometry and it 
is both interesting and legitimate to scrutinize the same situation in the case of large-scale 
electromagnetic fields. The purpose of the present study is to relax and complement the con- 
ventional approaches based on vacuum initial conditions where the length of the inflationary 
phase is assumed to be, a priori, largely immaterial for the physical and mathematical as- 
pects of the problem. The question we ought to address is, in short, the following: which 
are the initial conditions to be imposed on the evolution of the electromagnetic fields when, 
prior to the onset of inflation, a globally neutral plasma was present? In the case of the 
scalar and tensor modes of the geometry it suffices to contemplate, for some applications, 
the possibility of a non-vanishing number of phonons (or gravitons) in the initial state. In 
the electromagnetic case the situation is physically different since the electric and magnetic 
fields do not interact in the same way with the plasma. 

Consider a conformally flat Friedmann-Robertson- Walker (FRW) metric g^ u character- 
ized by the scale factor a(t). The evolution equations of the Hubble rate H = a/a are 

■2 



H 2 Ml = 



87T 

t 



Ptot + y + V(<p) 



HMp = -4tt 



(p 2 + (p tot + P tot ) 



(1) 



where ip denotes the inflaton field with potential V(<p) and the overdot denotes a derivation 
with respect to the cosmic time coordinate t; p tot and P tot are the covariantly conserved 
energy density and pressure of the plasma (i.e. d t p tot + SH (p tot + P tot ) = 0). The conformally 
flat metric g^ u can be explicitly written as g^uij) = a 2 (T)r]^ u where rj^ is the Minkowski 
metric and the conformal time coordinate r is related to the cosmic time t as a(r)dr = dt. 
In the case p tot = 3P to t we shall fix the attention on solutions of Eq. (1) with the following 
property 

where the parameter 5 < 1 measures the relative weight of the radiation background during 
the protoinflationary stage when the slow-roll parameter e = —H/H 2 is just about to drop 
below 1 and the second time derivative of the scale factor (i.e. a) is turning from negative 
to positive. In Eq. (1) and hereunder the physical variables (as opposed to their co moving 
counterparts) are denoted by a bar. The preinflationary (and hence decelerated) stage 
of expansion occurs for t <C H{ 1 . The protoinflationary regime corresponds instead to 
t ~ 0(Hi *); finally, for t 3> iff 1 the Universe inflates. 
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The total energy and pressure appearing in Eq. (1) depend on the charged and neutral 
species. The contribution of the charged species to the transport coefficients is computed 
under the hypothesis that the collisions between the particles of the same charge can be ne- 
glected (i.e. r± ^> T + and T± ^> as it happens for Lorentzian plasmas [3]. Furthermore, 
at high temperatures T± > H for interactions mediated by massless gauge bosons. Denoting 
with m± the masses of either the positive or negative charge carriers, their corresponding 
temperatures are the same as long as T± 3> m±. When the species become non-relativistic 
(i.e. T± < m±) the evolution of the temperature depends on the relative hierarchy of the 
charged and neutral concentrations. In the latter case the entropy conservation and the first 
principle of thermodynamics imply: 

— — — 47T 4 n 

d[a 2 {T + + T„)]+ 1 ad{aT r ) = Q ) 7 = — — ^, n + = n_ = n , (3) 

where, following the conventions of Eq. (1), n± = a 3 n± and n r = a 3 n r denote the comoving 
concentrations. During the preinflationary and the protoinflationary stages of expansion 
(see Eq. (2)) the radiation background dominates p tot (i.e. n r ^> n and 7 ^> 1). Equation 
(3) can be solved, to lowest order, for T + ~ T_ ~ T r and the three temperatures scale, 
approximately, as T ~ a~ K where re = (4 + 7)/(2 + 7) ~ 1 + 2/7 + 0(l/7 2 ). Note that 
7 oc ^asma where (7pi asm a = 1/(tio^d) ^ 1 is the plasma parameter quantifying, by definition, 
the inverse of the number of particles of charge q present in the Debye sphere (i.e. the sphere 
whose radius is given by the Debye length Ad = \JT / (87rn g 2 ) ). Provided g p iasma <^ 1 (i.e. 
when the Debye sphere contains a large number of particles) 7 is also much larger than 1; 
in a related context this requirement goes under the name of plasma approximation [4]. 

The contribution of the charged species to the Maxwell equations remains unsuppressed 
in the protoinflationary phase and, depending on the total number of efolds, during the early 
stages of the inflationary epoch. Indeed, for T ^> m± (i.e. T 3> m±a) the whole evolution 
of the charged species is Weyl invariant and this symmetry prevents the suppression of the 
conductivity in the relativistic regime. The simplest way to demonstrate the Weyl invariance 
of the full system rests on the explicit expression of the Vlasov-Landau equations for the 
distribution functions of the charged species written in a conformally flat metric of FRW 
type (see, e.g. [3]): 

df± 



Or 



+ v± ■ Vsf± =F q 



E + v±x B 



where E = a 2 S, B = a 2 B and v± = p/^Jp 2 + m±a 2 ] p denotes the comoving three- 
momentum. When the mass contribution can be neglected in comparison with p 2 Eq. (4) 
has the same form it would have in Minkowski space-time [3]. In the opposite situation (i.e. 
p 2 m\a 2 ) the mass term breaks explicitly Weyl invariance. The comoving concentrations 
of the charged species can be obtained by integrating the distribution functions over the 
comoving three-momenta so that the Maxwell equations are: 

V • E = 47rg(n+ V • B = 0, (5) 
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dB dE 
V x E = - — , V x B = 4irq(n + v+ - n^v.) + — . (6) 

Eqs. (4), (5) and (6) are then Weyl invariant in the relativistic regime and in the plasma 
approximation which guarantees an approximate common temperature of the charged and 
neutral species (see Eq. (3)). The latter result is quite known 2 and can be directly obtained, 
by means of simple scaling considerations, from the generally covariant form of Maxwell 
equations, i.e. 3 : 

d^V^D = 4«V=gj v , d^V^r) = o. (7) 

In four space-time dimensions both \/—gf^ u and yj—gj^ are invariant under the Weyl 
rescaling of the form g^x) = a 2 (x)ri f j iU where x labels a generic space-time coordinate. 
Consequently the whole system of Eq. (7) is Weyl invariant if y/^gj v is separately Weyl 
invariant. Consider, in this respect, a relativistic plasma with Ohmic current j v = af^u^ 
where w(x) is the physical conductivity. Since \/—gf fJ,u is separately Weyl invariant, the 
expression \J—gj v is Weyl invariant provided the combination au^ is invariant. Under Weyl 
rescaling transforms as — > w M = u^/a^x) because of g^ v u^u v = 1. For this reason 
\J—gj u is Weyl invariant provided the conductivity scales as a(x) — > a(x) = a(x)a(x). But 
this happens exactly in a relativistic plasma, where the physical conductivity scales as the 
first power of the physical temperature (i.e. w(x) oc T(x)) and it coincides, in practice, 
with the comoving conductivity a(x) oc T(x) = T(x)a(x) which is approximately constant 
whenever T ~ aT 1 . 

The lack of Weyl invariance of Eqs. (4)-(6) when T < m± is reflected in the general 
form of the conductivity; to simplify the expressions a hierarchy in the masses of the charge 
carriers can be assumed (for instance m + > m_ = m) so that, for a Lorentzian plasma [3], 
the conductivity is: 

a(a, 7 ) = - T } a,1 \ n , hm T(a, 7) oc a - 2/(7+2) = constant. (8) 



In the limit 7 ^> 1 the comoving temperature T = Ta is approximately constant; different 
situations, corresponding to specific values of 7, can be studied but they are less significant 
for the illustrative purposes of the present analysis. In the limit T ^> ma Eq. (8) implies 
a ~ T/q 2 (as it happens in the case of an relativistic plasma); in the opposite limit, a ~ 
T / q 2 ^1 T / {ma) . By subtracting the evolution equations of the velocities according to the 
standard procedure [4], it is possible to obtain the equation for the total current appearing 



2 This conclusion can be found within different perspectives in the book and in the papers reported in [5]. 

3 Note that and are, respectively, the Maxwell field strength and its dual; g denotes, as usual, the 
determinant of g^; in terms of the physical fields £ l and B % we have f l ° = £' l /a 2 and / u = — t % ^ k Bk/a 2 . 
Furthermore, as already mentioned after Eq. (4), comoving and physical fields are related as E = a 2 £ and 
B = a 2 B. 
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in Eq. (6): 



— + a(H + r±)J = -^(E + vx B) J = q(n + v + -n-v-), (9) 
Or An 

where terms like JxB have been neglected and where v = (m + v + +m-V-) / (m + +m_). Since 
r± > H the total Ohmic current becomes indeed J ~ a(E + v x B) with a = u^/ (47rar±). 

Equation (8) can be phrased in terms of the number of inflationary efolds iV = In (a/ai) 
where a\ denotes the scale factor at the onset of the inflationary phase; thus we have a(N) = 
<7*/y/i + e N ~ N *. For N < N* the conductivity is approximately constant while for iV > N* 
it is exponentially suppressed. To leading order in g p i asma the critical number of efolds 
A* depends on the temperature reached during the protoinflationary phase and it can be 
estimated by recalling that at the onset of the inflationary phase S < 1 in Eq. (2). Requiring 
that, at most, 5 ~ 1 an upper bound on N* can be obtained: 

1 1 / m \ 

N* = -0.253 + - In £ - - In g th - In ^— J , (10) 

where g t h denotes the effective number of relativistic degrees of freedom; £ = Hi/M P = 
y/neAn where An = 2.43 x 10~ 9 is the amplitude of the scalar power spectrum at the pivot 
scale k p = 0.002 Mpc -1 and e is the slow-roll parameter introduced after Eq. (2). For the 
numerical estimates we shall adopt the values of the cosmological parameters determined in 
terms of the WMAP data alone in the minimal concordance model (see last three papers 
of Ref. [6] for the latest release). Since the maximal protoinflationary temperature should 
not exceed the energy density of the background, we have that, at most, iV* ~ 36.78 — 
0.25 In (^th/100) + 0.5 In (£/10~ 5 ) - In (m/GeV). The value of iV* can be compared with the 
maximal number of efolds presently accessible by large-scale observations, i.e. 

Nl = 62 .2 + i In (*) - In (!*) + I In ( ^ ,V (11) 
2 VIO-V V0.7/ 4 V4.15xl0- 5 /' v ; 

which is close, by construction, to the minimal number of efolds iV m i n needed to solve the 
kinematic problems of the standard cosmological model (i.e. iV min ~ Aq). Recalling the 
fiducial set of cosmological parameters quoted before [6], Eq. (11) gives Aq = 63.6 + 0.25 In e. 
Equation (11) is derived in the sudden reheating approximation but Aq can be larger if right 
after inflation the Universe expands at a rate which is slower than radiation down to the 
big-bang nucleosynthesis curvature scale. In the latter case the estimate for Aq increases by 
(9(14) efolds so that Aq ~ 78.3 + 0.33 In e in agreement with previous estimates [7]. 

If A^ ~ Atot ^> A 7 * + Amin, Weyl invariance is broken before the onset of the last Aq ~ 
(9(63) efolds of inflationary expansion; the electromagnetic fields are normalized when the the 
protoinflationary conductivity is suppressed as e -( N t°t- N *)/ 2 . But this means that the sources 
do not contribute, in practice, to the initial conditions which are accurately fixed by quantum 
mechanics. Conversely if A" min < N < N* + N miD the conductivity will be constant for the 
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first iV* efolds and then it will be exponentially suppressed as e { N *- N ^)/ 2 (if iV ~ N min ) and 
as e~ Nmin ^ 2 (if N ~ iV miTi + iV*). While the point here is not to endorse a specific duration 
of the inflationary phase, when A^" min < N < N* + N min the last (9(63) efolds of inflationary 
expansion may start when the conductivity did not undergo a substantial suppression. For 
instance, when the mass range of the lightest charge carrier is O(GeV) and if N ~ 0(N min ) 
the conductivity is still almost constant (9(30) efolds prior to the end of inflation. In this class 
of physical situations the normalization of the electric and magnetic fields does not follow 
from the quantum mechanical initial conditions but rather from the standard conducting 
initial conditions [4]. From Eqs. (5)-(6) and (9) E is solenoidal (because of the global 
charge neutrality) but also J must be solenoidal since in the plasma rest frame E ~ J /a. 
In the regime of high conductivity the displacement current can be neglected and therefore 
the appropriate initial conditions for the electromagnetic fields at r* are simply given by 

B(x,t*) = B^\x), E(x,t.) = — a ^, (12) 

47T a * 

— > — * — > — * — > — * 

together with the conditions V-E — V-B — V-J — 0. Let us now characterize, for 
immediate convenience, the electric and magnetic fields by means of their associated power 
spectra defined, in Fourier space, as 

2tt 2 

(BiikrfBj&r)) = -^P B (Kr)P ij (k)5^(k+^, (13) 

where Pij(k) = 5ij — kikj is the transverse projector; exactly the same definition holds in the 
case of the electric fields whose related power spectrum will be denoted by Pe(&, t). 

The vacuum and the conducting initial conditions lead to different power spectra which 
can be compared. Consider, for instance, the class of models where the gauge coupling is 
dynamical during inflation according to the action / d 4 x i/— 9 VA / '^f^ u / (167r) [8]. When 
A increases the gauge coupling decreases and this possibility fits with the presence of a 
curvature singularity in the past history of the inflationary dynamics: since the gravitational 
coupling gets strong, it is also appropriate to contemplate the case when a/A is initially 0(1), 
it increases during the inflationary phase as ~ (— t) x I 2 ~ v and it decreases again during 
the reheating process (see last four articles of Ref. [8]). The Hamiltonian describing the 
evolution of the classical electromagnetic inhomogeneities is given, in Fourier space, by 



(14) 



where the prime denotes a derivation with respect to the conformal time coordinate r. In 
the Coulomb gauge the (comoving) electromagnetic fields B = V x y and E = —n; Eq. (14) 
is invariant under the duality transformation \/\ — > -4= and — > —ky^, y^ — > |7T£ (see 
first two references in Ref. [8]). The Hamilton equations derived from Eq. (14) become: 
yr — ^k+ [(V^)'/V^]j7jfc> an d = ~^ 2 Vk ~~ [(v^y/v^A]^;- In the case of the amplification 
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of vacuum fluctuations the power spectra of the magnetic field have a spectral slope /c nB_1 
with fie = (6 — 2z/); in scale invariant limit (i.e. v = 5/2) the present value of the power 
spectrum can be estimated as [8] 



log [JP B (k, r )/Gauss] = -10.85 + 0.5 



log 



A 



n 



2.43 x 10" 9 



+ log 



RO 



4.15 x 10" 5 



(15) 



where P B (k,r) = a 4 (r)Pg(A;, r). The result of Eq. (15) assumes, for consistency with Eqs. 
(10) and (11), the sudden reheating approximation. When the modes of the field are inside 
the Hubble radius during inflation (i.e. kr > 1) the electric and the magnetic power spectra 
are of the same order (i.e. P E {k,T) ~ P B (k,T)). In the opposite case (i.e. kr <C 1) the 
electric power spectra are suppressed, in comparison with their magnetic counterpart, as 
P E (k,r) ~ \kr\ 2 P B {k, r). In the sudden reheating approximation the conductivity raises 
again at the end of inflation and this step further suppresses exponentially the electric fields. 

The estimate of Eq. (15) can be compared with magnetic field obtainable in the frame- 
work of the same toy model but in the case of the conducting initial conditions. For sake 
of concreteness we will have that for r < — r* the conductivity and A will both be constant 
(7 = 0"* and y/X = v^*- For r > — t* the conductivity will vanish, d T y/\ > and we can 
assume, for sake of comparison with the vacuum case of Eq. (15), that \f~\ ~ (— t) 1 ^ 2 ^". 
Therefore, the electromagnetic fields subjected to the conducting initial conditions (12) and 
subsequently amplifield by the evolution of v^A are: 



E(k, t, r*) = —ik x B(k, r*) Q(x, z), B(k, r, r*) = B(k, r*) T(x, z), 



(16) 



where x = kr* and z = —kr; recalling the definitions of Hankel functions of first and 
second kind (i.e. and H^) it is useful to define the following combination Pj) l \z) = 

2vHW(z) - zH^z) (and similarly for PP(z) = 2uH^(z) - zH^zj). Consequently the 
functions Q(x,z) and T{x,z) appearing in Eq. (16) can be written as: 



G(x,z) = 



F(x,z) = 



m 
i 

V 
in 



if 1 

Z I X 



Pl 1 \z)Pl 2 \x) - Pl 1 \x)Pl 2 \z) 



H^\x)P^\z)-H^\x)Pj ) 2 \; 



2v — TjX 



Hj 2 \x)Hi 1 \z) - Hi 1 \x)Hl 2 \z) 



+ X 



Hl%(x)HW(z) 



(17) 



where 77 = k/Ancr*. By setting r = — r* the expressions of Eq. (17) reproduce the conducting 
initial conditions given in Eq. (12). Furthermore, the associated power spectra can be writ- 
ten, respectively, as P B (k,T,T*) = P B (k,T*) IJ 7 ^, x)| 2 , and P E (k,T,T*) = P B (k,T*) \Q(z, x)\ 2 . 
In the limit |r*| |r e | (and x < 1, z < 1) the electric power spectra are always exponen- 
tially suppressed already during inflation in comparison with the magnetic power spectra. 
From Eq. (17), at the end of the inflationary phase, P B (k, r e , r*) ~ Psik, r*)(a e /a*) 2u ~ 1 and 
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P E (k,r e ,n) ~ r)P B (k,n) (a e /a*y- 2u . Since rj = A;/(47ra) ~ fc/T, < 1 and 1/2 < v < 5/2 
the suppression of the electric fields is always much larger than in the case of vacuum initial 
conditions. To assess the amplitude of the magnetic fields we shall consider, for sake of 
comparison with Eq. (15) the case v — 5/2. To avoid supplementary assumptions, the am- 
plitude of the magnetic power spectrum during the protoinflationary phase can be set to its 
maximum value compatible with the closure bounds. In the latter case the power spectrum 
of the magnetic field will be given by 

log[^P B (fc,T,,T )/Gauss] = -60.15 + 0.5 log + (2a - 1) log (J^) 

+ 0.43[(^ — 1/2)(JV — JV*) — 2(N — Ni)], (18) 

where A = ^JPs(k, T\)/ (HfMp) < 1 measures the fraction of energy density stored in the 
magnetic field at T\\ H r accounts for the possibility of a delayed radiation-dominated phase 
between the end of inflation and the onset of big-bang nucleosynthesis. The exponent a 
depends on the expansion rate between the end of the inflationary phase and the onset of 
the standard (i.e. post-inflationary) radiation-dominated epoch. Equation (18) has several 
interesting limits. Recalling that the the interval of v is restricted, from kinematical consid- 
erations, tol/2<z/<5/2, in the case v — 1/2 there is no amplification due to the evolution 
of the gauge coupling and therefore the upper bound on the magnetic field intensity is around 
(9(10~ 61 ) Gauss in the sudden reheating approximation where H r ~ H e . This is simply the 
magnetic field one would obtain from the protoinflationary initial conditions assuming the 
minimal amount of inflationary efolds. Conversely, if v — 5/2 and iV = N 1 ~ (9(65) and 
N* ~ (9(35) the maximal magnetic field turns out to be 10~ 35 Gauss (in the sudden reheating 
approximation) which can become of the order of O(l0~ 23 ) Gauss for a stiff post-inflationary 
phase extending down to the nucleosynthesis scale. 

The implications of a globally neutral plasma during the protoinflationary stage of expan- 
sion have been investigated. This idea did not receive specific attention so far even if it seems 
rather natural in the light of similar attempts carried on in the case of the conventional scalar 
and tensor modes of the geometry. If the plasma is relativistic, Weyl invariance prevents 
the suppression of the conductivity which starts being diluted as soon as Weyl invariance 
gets broken by the masses of the charge carriers. The presence of a relativistic plasma in 
the initial conditions of the inflationary dynamics changes qualitatively and quantitatively 
the initial data to be imposed for the evolution of the large-scale electromagnetic inhomo- 
geneities. The conducting initial conditions for the amplification of large-scale magnetic 
fields have been contrasted with the conventional vacuum initial conditions. In a class of 
specific examples the amplification of the magnetic fields and the suppression of the electric 
fields has been shown to depend explicitly on the number of efolds of the inflationary phase. 
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